UNIVERSALITY OF KPZ EQUATION 
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Abstract. We introduce the notion of energy solutions of the KPZ equation. 
Under minimal assumptions, we prove that the density fluctuations of one- 
dimensional, weakly asymmetric, conservative particle systems with respect 
to the stationary states are given by energy solutions of the KPZ equation. As 
a consequence, we prove that the Cole-Hofp solutions are also energy solutions 
of the KPZ equation. 



1. Introduction 

In a seminal paper 11331 . Kardar, Parisi and Zhang proposed a phenomeno- 
logical model for the stochastic evolution of the profile of a growing interface 
ht{x). The so-called KPZ equation has the following form in one dimension: 

dth = DA/i + a{Vhf + aWt, 

where Wt is a space-time white noise and the constants D, a, a are related 
to some thermodynamic properties of the interface. The quantity ht (x) rep- 
resents the height of the interface at the point x e M. From a mathematical 
point of view, this equation is ill-posed, since the solutions are expected to look 
locally like a Brownian motion, and in this case the nonlinear term does not 
make sense, at least not in a classical sense. This equation can be solved at a 
formal level using the Cole-Hopf transformation Zt{x) = cxp{a/Dht{x)}, which 
transforms this equation into the stochastic heat equation 

dtZ = DAZ + acj/DZWt- 

This equation is now linear, and mild solutions can be easily constructed. We 
will call these solutions Cole-Hopf solutions. It is widely believed that the 
physically relevant solutions of the KPZ equation are the Cole-Hopf solutions. 
However, the KPZ equation has been so resistant to any attempt to mathemat- 
ical rigor, that up to now it has not even been proved that Cole-Hopf solutions 
satisfy the KPZ equation in any meaningful sense. Some interpretations that 
allow rigorous results are proved to give non-physical solutions [13|. Up to our 
knowledge, the best effort in this direction corresponds to the work of Bertini 
and Giacomin [7|. In that work, the authors prove two results. First they prove 
that the Cole-Hopf solutions can be obtained as the limit of a sequence of mol- 
lified versions of the KPZ equation. And secondly, they proved that the Cole- 
Hopf solution appears as the scaling limit of the fluctuations of the current for 
the weakly asymmetric simple exclusion process (WASEP), giving mathemati- 
cal support to the physical relevance of the Cole-Hopf solution. 
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In dimension d = I, a conservative version of the KPZ equation can be ob- 
tained by defining yt ^Vht: 

This equation has (always at a formal level!) a spatial white noise as an in- 
variant solution. In this case is even clearer that some procedure is needed in 
order to define yf in a proper way. 

Since the groundbreaking works of Johansson fSOl and Baik, Deift and Jo- 
hansson 151, a new approach to the analysis of the KPZ equation has emerged. 
The general strategy is to describe various functionals of one-dimensional asym- 
metric, conservative systems in terms of determinantal formulas. These deter- 
minantal formulas turn out to be related to different scaling limits appearing 
in random matrix theory. We refer to the expository work 112 111 for further ref- 
erences and more detailed comments, and to [42|, [1| for recent advances. 

Another approach to the analysis of fluctuations of one-dimensional conser- 
vative systems was proposed by Balazs and Seppalainen |[9|. They call their 
approach microscopic concavity I convexity and it is exploited in fHl in order to 
prove that the Cole-Hopf solution of the KPZ equation has the scaling expo- 
nents predicted by physicists. 

The main drawback of all these approaches is the lack of robustness. The 
microscopic Cole-Hopf transformation used in | 7 | works only for the WASEP, 
and any other interaction different from the exclusion principle rules out this 
approach. Determinantal processes appear in a natural way for particle sys- 
tems that can be described by non-intersecting paths. The microscopic concav- 
ity/convexity property can be defined only for attractive systems, and it has 
been proved only for the asymmetric simple exclusion process (ASEP) and the 
totally asymmetric, nearest-neighbors zero-range process under very restric- 
tive assumptions. Moreover, up to our knowledge all the approaches to the 
KPZ equation go through Bertini and Giacomin construction of the Cole-Hopf 
solutions (see however 139]). 

It is widely believed in the physics community that the KPZ equation gov- 
erns the large-scale properties of one-dimensional, asymmetric, conservative 
systems in great generality. The microscopic details of each model should only 
appear through the values of the constants D, a and a. In this article we pro- 
vide a new approach which is robust enough to apply for a wide family of one- 
dimensional, asymmetric systems. The payback of such a general approach 
comes at the level of the results: we are not able to prove the precise results of 
Bertini and Giacomin, and we can not recover the detailed results obtained by 
the random matrix theory approach. However, our approach is robust enough 
to give information about sample path properties of the solutions of the KPZ 
equation. 

As a stochastic partial differential equation, the main problem with the KPZ 
equation is the definition of the square y^. Spatial white noise is an invariant 
solution of the equation, and it is expected that physically relevant solutions 
look locally like white noise. 

Our first contribution is the notion of energy solutions of the KPZ equation 
(see Section |2^ . Various attempts to rigorously define a solution of the KPZ 
equation have been made. One possibility is to regularize the noise Wt and 
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then to turn the regularization off. For the regularized problem, several prop- 
erties, like well-posedness of the Cauchy problem and existence of invariant 
measures can be proved [4r|, |181. However, the available results hold in a 
window which is still far from the white noise Wt. Another possibility is to 
regularize the nonlinearity Vyf II16L One more time, this procedure gives 
well-posedness in a window which is far from the KPZ equation. Yet another 
possibility corresponds to define the nonlinear term through a sort of Wick 
renormalization ||28]| . lH. However, this procedure does not lead to solutions 
with the right scaling properties \ TS\. Our notion of energy solutions is strong 
enough to imply some regularity properties of the solutions which allow to jus- 
tify some formal manipulations. 

We introduce the notion of energy solutions of the KPZ equation in order 
to state in a rigorous way our second contribution. Take a one-dimensional, 
weakly asymmetric conservative particle system and consider the rescaled 
space-time fluctuations of the density field 3^" (see the definition in Section 
12.21 1. The strength of the asymmetry is of order 1/y/n. For the speed-change 
simple exclusion process considered in II19II and starting from an stationary dis- 
tribution, we prove that any limit point of 3^" is an energy solution of the KPZ 
equation. The only ingredients needed in order to prove this result are a sharp 
estimate on the spectral gap of the dynamics of the particle system restricted 
to finite boxes (stated in Proposition 13. 4D and a strong form of the equivalence 
of ensembles for the stationary distribution (see Proposition 13. 7| >. Therefore, 
our approach works, modulo technical modifications, for any one-dimensional, 
weakly asymmetric conservative particle system satisfying these two proper- 
ties. In particular, our approach is suitable to treat models like the zero-range 
process and Ginzburg-Landau model in dimension d = 1, for which the meth- 
ods mentioned before fail dramatically. Our approach also works for models 
with finite-range, non-nearest neighbor interactions with basically notational 
modifications. 

We consider in this article speed-change exclusion processes satisfying the 
so-called ^radie^i^ condition. Notice that we need to know the invariant mea- 
sures of the model in order to state the equivalence of ensembles. It has been 
proved Ii37| for the speed-change exclusion process that the invariance of a 
Gibbs measure under the symmetric dynamics is preserved by introducing 
an asymmetry, if and only if the model satisfies the gradient condition. It 
is only at this point that we need the gradient condition. In particular, our 
approach also works for weakly asymmetric systems for which the invariant 
measures are known explicitly, even if the gradient condition is not satisfied. 
In view of this discussion, we say that energy solutions of the KPZ equation 
are universal, in the sense that they arise as the scaling limit of the density 
in one-dimensional, weakly asymmetric conservative systems satisfying fairly 
general, minimal assumptions. 

In order to prove this theorem, we introduce a new mathematical tool, which 
we call second-order Boltzmann-Gibbs principle. The usual Boltzmann-Gibbs 
principle, introduced in [10| and proved in [W\ in our context, basically states 
that the space-time fluctuations of any field associated to a conservative model 
can be written as a linear functional of the density field y^. Our second-order 
Boltzmann-Gibbs principle states that the first-order correction of this limit 
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is given by a singular, quadratic functional of the density field. It has been 
proved that in dimension d > 3, this first order correction is given by a white 
noise [14]. It is conjectured that this is also the case in dimension d ~ 2 and in 
dimension d = 1 if our first-order correction is null. 

The rest of this paper is organized as follows. In Section |2] we give precise 
definitions of the model considered here and we state the results proved in the 
rest of the article. In Section |2?T] we define the speed-change exclusion process 
and we state some of its basic properties. In Section [2l2] we give an overview of 
various scaling limits of the density of particles for the speed-change, simple 
exclusion process, and in particular we state Bertini and Giacomin's result. 
In Section 12.31 we give rigorous definitions of what we understand by weak 
and energy solutions of the KPZ equation and we state our main result. In 
few words, our main result states that the density fluctuation field is tight, 
and any limit point is an energy solution of the KPZ equation. We state two 
corollaries, the first one gives the Holder exponent of the energy solutions of 
the KPZ equation. The second gives an answer of the following open problem in 
Q: does the Cole-Hopf solution of the KPZ equation actually satisfies the KPZ 
equation in any meaningful way? Applying our main result to the WASEP, we 
prove that the stationary Cole-Hopf solution is an energy solution of the KPZ 
equation. In Section 12.41 we define the stochastic growth model associated to 
the speed-change, simple exclusion process and we restate the definitions of 
Section in terms of this growth model. We prove a central limit theorem 
for the current through a bond. For simplicity, we assume that the drift has 
zero average, but our results remain true if we look at the current across a 
characteristic. Starting from this result we prove a convergence result for the 
height fluctuations on this stochastic growth model. 

In Section [3] we state and prove the second-order Boltzmann-Gibbs prin- 
ciple, which is is the main technical innovation of this article. In Section 
13.11 we review Kipnis-Varadhan and spectral gap inequalities and the equiv- 
alence of ensembles, and we state Propositions 13.51 and 13.71 We point out 
here that the equivalence of ensembles we need is of second-order and in con- 
sequence it is finer than the result usually found in the literature. For the 
reader's convenience we give in Appendix|A|the proof of this result in our sim- 
ple case of Bernoulli uniform measures. In Section 13.21 we prove the second- 
order Boltzmann-Gibbs principle, only relying in Propositions l3.5l and l3.71 The 
proof follows from a multiscale analysis introduced in [26|. The multiscale 
analysis has two steps. The first one, which we call the seed is reminiscent of 
the well-known one-block estimate. In the second step we apply a key iterative 
bound to go from a microscopically big block to a macroscopically small block. 

In Section|4]we prove Theorem l2.4l The proof follows the classical scheme to 
prove convergence theorems in probability. In Section 143] we prove tightness 
of the density fluctuation fields and in Section 14.21 we prove that any limit 
point is a stationary energy solution of the KPZ equation. We conjecture that 
energy solutionst starting from the stationary state are unique in distribution. 
Conditioned to this uniqueness result, convergence follows. 

In Section [5] we prove the convergence results for the current and height 
fluctuation fields. The central limit theorem for the current follows from an 
idea of Rost and Vares 1 401 and we follow the approach of [i29il and ||26J . The 
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height fluctuation field formally corresponds to the integral of the density fluc- 
tuation field. The convergence of the height fluctuation field does not follow 
directly from the convergence of the density, since we need to deal with the 
constant of integration, which is a non-trivial process that we relate with a 
sort of mollified current process. 

2. Notation and results 

2.1. The model. Let fl = {0,1}^ be the state space of a continuous-time 
Markov chain r/t which we will define as follows. We say that a function 
/ : -> R is local if there exists R = R{f) > such that f{ri) = f{() for 
any r;,^ e such that r]{x) = ^{x) whenever > R. In other words, we say 
that / is local if f{r]) depends only on a finite number of coordinates of r] (in 
this case at most 2i?+ 1). Let c : O — > M be a non-negative function. We assume 
the following conditions on c: 

i) Ellipticity: There exists eo > such that eo < c{ri) < e^^ for any 77 e fi. 

ii) Finite range: The function c(-) is local. 

iii) Reversibility: For any 77, ^ G such that r]{x) = S,{x) whenever x ^ 0,1, 
c{v) = c(0. 

For any x £ Z\et : fl ^ Qhe the translation in x: Txri{z) = ri{x + z) for 
any rj £ Q and any z e Z. For a function / : — > R we define r^f : fl ^ M as 
Txfiil) — f{'''x'ri) foi" any j] e Q. We will also assume a fourth condition, which is 
the most restrictive one: 

iv) Gradient condition: There exists a local function h : fl ^ fl such that 

c{t]){t]{1) - 77(0)) = Tih{T]) - h{7]) for any tj e Q. 

The exclusion process (possibly asymmetric) with speed change is defined 
as the Markov process {77"; t > 0} generated by the operator L„, whose action 
over local functions / : — > R is given by 

Lnfiv) = "-^ X! ^x{v){PnV{x){l - + 1)) + QnVix + 1)(1 - Vi^))}'^ x.x + 1 f iv) , 

xez 

where n € nQ, c^{ri) = T^c{ri), V^.x+ifiv) = ''^^^) - fiv), Pn and g„ are non- 
negative constants such that p„ + 9„ = 1 (eventually we will choose p„ — Qn ~ 
a/ yjn with a 7^ 0) and jf^^^^ is given by 

{77(x + 1), z = X 
r]{x), z^x + 1 
ri{z), z ^ x,x + 1. 

In order to explain the meaning of conditions i)-iv) let us assume by now 
that Pn — q-n — 1/2. In this case the process 77" is said to be symmetric. Con- 
dition i) ensures that the process is well defined (see Chapter 1 of 135] for a 
comprehensive discussion) for any choice of p„, g„. For p £ [0, 1] let Vp be the 
Bernoulli product measure in of parameter p. This means that for any two 
finite, disjoint sets A,B £ Z, 

i^M^) - 1 Va; e A, 7/(77) = V77 e B) = pl^l (1 - 



Here and below we use the convention N = {1,2,...} 
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where \A\, \B\ denote the cardinahty of the sets A, B respectively. Under condi- 
tion iii), the measures {vp-. p £ [0, 1]} are invariant and reversible with respect 
to the evolution of r/". Under condition i), these measures are also ergodic with 
respect to the evolution of r/". Condition i) also tells us that the dynamics of 77" 
is comparable to the dynamics of the simple exclusion process without speed 
change: let Lff be the generator associated to the choice c( ) = 1. Then, 

j f{-Lr,f)dyp < J f{-Ll-f)dyp < e^'J f(-L^f)dvp (2.1) 

for any local function /. We will return to the meaning of this bound later. 

When the process 77" is asymmetric (that is, when p„ ^ qn), it is not true in 
general that the measures Up are invariant with respect to the evolution of 
In fact, according to ||37]| . the family {vp;p G [0, 1]} is invariant with respect 
to the evolution of ?/" if and only if the condition iv) is satisfied. It is exactly 
for this reason that we assume the restrictive condition iv). In this case, the 
measures Vp are no longer reversible with respect to the dynamics. 

2.2. Scaling limits. In this section we recall various scaling limits previously 
obtained for the density of particles with respect to the process 77". These 
results are known in the literature as hydrodynamic limits. First we recall 
a law of large numbers for the density of particles. Let ■n^{dx) be the positive 
measure in M defined as 

Kidx) = i V?7("(a;)(5^/„(da;), 
n ^ — ' 

where (5^/„(dx) is the Dirac mass at x/n e R. The process 7r"(c?2:) is known 
as the empirical density measure associate to the process 77". When the dis- 
tribution of 77q is equal to Vp with p e [0, 1] then the distribution of 77" is also 
equal to Vp for any later time i > 0. In particular, the sequence of measures 
{■K^{dx)]n e N} converges in probability to the measure pdx when 71 — > 00 for 
any t > 0. 

Let uq : R ^ [0, 1] be a strictly positive, continuous by parts function such 
that there exists p € (0, 1) for which / |Mo(a;) — p\dx < +00. Let {^"; n e N} be a 
sequence of Bernoulli product measures in fl, defined by the relation p^^{r]{x) — 
1) = UQ{x/n). Let us denote by P^^. the distribution of the process 77" with initial 
distribution /i". The following theorem is known as the hydrodynamic limit of 
the process 77": 

Proposition 2.1 ( Il27l [32] ). Let us take pn — qn = a/n with a e R and n big 
enough. Let uq :R ^ [0, 1] and {^"; n £N} be as above. Then n^{dx) converges 
in probability to the measure u{t,x)dx with respect to P^™, where {u{t,x);t > 
0, a; e M} is the solution of the hydrodynamic equation 



(2.2) 



\dtu = l/2A^(u) -aV/3(M) 
[-u(0, •) =uo{x) 

and ip{p) = J hdvp, = xip) J cdvp. 

The Einstein relation, which holds in great generality, states that Pip) = 
x{p)^'{p). The quantity x{p) — xif) is known in the literature as the conduc- 
tivity or susceptibility of the system. As we said before, the gradient condition 
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iv) is very restrictive. From the point of view of modeling, this condition is not 
so restrictive. In fact, for any m e N there exists a choice for 0(77) such that 
(p{p) ~ p™ (see 1 25 1). By hnearity, for any polynomial q{p) there are a constant 
K and rate 0(77) such that (p{p) ~ q{p) + K. 

The symmetric case a = in Proposition [211] was treated in II27II for a model 
with interactions of Ginzburg-Landau type, and the weakly asymmetric case 
was treated in |32|. In both cases, the gradient condition iv) is fundamental. 
The non-gradient method developed in Il38l |43]| allows to generalize this theo- 
rem (in the case a = 0) to the case on which iv) is not satisfied. In that case the 
function ip{p) is given in terms of a variational formula, and in particular ip{p) 
is not explicit. Under the additional hypothesis <^( ) e C^((0, 1)), Il23]| proved 
Proposition 12.11 when a = for the very same model considered here. Later, 
in [6] it is proved that Lp{-) e C°°((0, 1)), closing the proof of Proposition 12 . ll for 
this model. 

In this article we are interested in a central limit theorem for the density 
of particles, with respect to the hydrodynamic limit considered above. The 
current state of the art restrict ourselves to the equilibrium situation, that is, 
when the initial distribution of 77" is equal to Vp for some p € (0, 1) (the cases 
p = 0, 1 being trivial). Let us fix now and for the rest of the paper a density 
p e (0, 1) and let P„ be the distribution of 77" with initial condition Vp. We 
denote by E„ the expectation with respect to P„. Let S{M) be the Schwartz 
space of test functions and let S' (K) be the space of tempered distributions 
in M, which corresponds to the topological dual of 5(K) (we give more precise 
definitions in Sect. [2^. The fluctuation field {yt °;t > 0} is defined as the 
5'(R)-valued process given by 

y-'^^iG) ^ j=Y,{ri\x) - p)Gix/n) 

for any function G e S{R). Calculating the characteristic function of the ran- 
dom variable X' ^CC), it is easy to see that for any fixed time t>0, the process 
y"''^ converges in distribution to a spatial white noise of variance x(p)- For the 
sequence of processes {3^" °; n G N} the scaling limit is the following: 

Proposition 2.2 (|jl5l [nl [22l), Let us take Pn — Qn = a/n with a E M. and n 
big enough. The sequence {y^'°;n € N} converges in distribution with respect 
to the Ji-Skorohod topology inV^O, 00), 5' (R)) to the process 3^", solution of the 
Ornstein-Uhlenbeck equation 

dy? = lv'{p)Ay?dt ~ af3'{p)\/y°dt + ^x{p)^'{p)^dWu (2.3) 
where Wt is a space-time white noise of unit variance. 

The proof of this result relies on a replacement known as the Botzmann- 
Gibbs principle, which was introduced by H. Rost lllOl . When the gradient 
condition iv) is satisfied, this theorem was proved by IfTSlI for the case a = 0. 
Later a more robust proof was given by IdlJ , which can be extended to the case 
a ^ 0. When the gradient condition iv) is not satisfied, a sketch of the proof is 
given in [I22II . following the method outlined in [12]. 

The drift term in ( |2.3| l is irrelevant in the following sense. Performing a 
Galilean coordinate transformation, the drift term in equation ( 12. 3D can be 
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removed. This is clear if /3'{p) is equal to 0, in which case the equation ( 12. 31 * 
does not depend on a. In other words, the drift term does not give rise to 
stochastic fluctuations of the density field. In a more precise way, let us define 
the modified fluctuation field 3^" as 

where v{p) = aj3'{p) is the velocity associated to the system. Starting from 
Proposition 12 .21 we see that the sequence n e N} converges to the process 
yt, solution of the equation 

dyt - ]^ip'{p)Aytdt + ^Jx{pW{p)^dWu (2.4) 

which corresponds to equation ( 12.31 1 with a = 0. 

As we mentioned before, the asymmetry is too weak in order to induce a 
stochastic fluctuation in the density of particles. According to Q, a non triv- 
ial fluctuation due to the asymmetry appears when p„ — q„ = a/y^, on which 
case the limiting process yt has a qualitatively different evolution from the 
Ornstein-Uhlenbeck process (solution of ( I2.4I >). For this choice of the asymme- 
try, the modified fluctuation field is given by 

XXG) = ^ E (''r(^) - p)G{x/n - v{p)tn^'% 

Proposition 2.3 (17)). Let us take pn — Qn = a/^/n and c = 1. Let uq : M. ^ 
[0, 1] be a continuous function and let /i" be the measure appearing in Proposi- 
tion \2.1\ Then the process 3^" converges in distribution with respect to P^n to 
the Cole-Hopf solution of the KPZ equation 

dyt = ^Aytdt + aWyfdt + ^ut{l-ut)\/dWt (2.5) 

with initial condition W^", where ut is the solution of the hydrodynamic equa- 
tion ( 12. 2D and W"° is a Gaussian process in M of mean zero and covariance given 

by uo{x){l - uo{x))S{x,y). 

From now on and up to the end of the paper, we take p„ — q„ = a/-/n and for 
ease of notation we assume that v{p) — 0. In Section |2^ we will explain better 
what do we call a Cole-Hopf solution of equation ( 12.51 1. The KPZ equation was 
introduced in |33 | as a continuum model for surface growth. The proof of this 
theorem makes explicit use of the fact that the simple exclusion process is a 
totally integrable system. In a more precise way, in lEl the authors exploit 
a non-linear transformation of the process, discovered by Gartner [24| which 
linearizes the evolution of 3^". This transformation is a microscopic analogous 
of the Cole-Hopf transformation, and it is of help only in the case c = 1. 

The main result of this article is a generalization of Proposition 12.31 in the 
equilibrium case. But in order to write the result in a precise way, we need to 
introduce various definitions. 

2.3. KPZ equation. The KPZ equation, which was introduced by Kardar, 
Parisi and Zhang in the celebrated paper 1331 formally reads 

dht = ^Ahtdt-^^^{Vht)^dt + ^x{p)v'ip)dWt, (2.6) 
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where ht is a stochastic process with values on the set C(R) of continuous func- 
tions in M and Wt is a space-time white noise, that is, a Gaussian process of 
mean zero and covariance S{x — x')6{t — t'). Formally defining yt = Vht we 
obtain the conservative KPZ equation: 

dyt - ^^ytdt - ^^^vyfdt + ^x{pW{p)ydWt. (2.7) 

In the literature this equation is sometimes called stochastic Burgers equa- 
tion. The stochastic Burgers equation has been studied in detail for the case 
on which the noise has the strength enough to regularize the solutions of the 
equation. Therefore, we prefer to reserve the term stochastic Burgers equa- 
tions for the well studied cases on which the noise regularizes the equation. 
It turns out that, at least formally, the spatial white noise of variance xip) is 
invariant under the evolution of 02. 7D . Let C^(M) be the space of infinitely dif- 
ferentiable functions / : R M of compact support. For each / , m e N we define 
in (R) the norms 

||/||(,,„)=sup|x'/(")(x)|, 

where /(™) is the m-th derivative of /. Now we define the distance ds : C^(R) x 
C^{R) ^Ras 

;,meN 

for any f,g e C^(R). The Schwartz space S{R) is defined as the closure of 
C^(R) with respect to the distance ds{-,-)- The space 5(R) coincides with the 
set of infinitely differentiable functions / : R ^ R such that ||/|| (;.,„) < +oo 
for any Z,to e N. The space of tempered distributions S'(R) is defined as the 
topological dual of 5(R). For T e (0,oo), let C([0,T],5'(R)) be the space of 
continuous functions from [0, T] to 5'(R). 

The space C([0,T],5'(R)) is the space on which the solutions of the KPZ 
equation i2.7i will live. For e > we define : R M by 

ie{x){y) = e"H(.T <y <x + e). 

We say that a process {yut e [0,r]} with trajectories in C([0,T],5'(R)) and 
adapted to some natural filtration {Tt; t e [0, T]} is a weak solution of the KPZ 
equation ( |2.7| t if: 

i) There exists a process {At;te [0, T]} with trajectories in C([0, T], 5'(R)) 
and adapted to { J"*; t e [0, T]} such that for any G e S{R), 

lim r / y,{t,{x)f^^^^^±^l^-^dxds^At{G). (2.8) 

ii) For any function G e 5(R) the process 

Mt{G) = yt{G) - yo{G) - ^ ys{G")ds - ^^^mg) (2.9) 

is a martingale of quadratic variation x{p)y^' {p)t\\G\\i- 
Here and below we write ||G||f = / G'{x)^dx. This notion of solution of the 
KPZ equation is new. Moreover, it seems to us that any effort to define a mean- 
ingful notion of solution of ( I2.7D has not been well succeeded. Now we will 
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introduce a stronger notion of solution, which captures well some of the partic- 
ularities of the solutions of i2.7i . Let {3^t; t e [0, T]} be a weak solution of ( I2.7I ). 
For < s < t < T, let us define the fields 

IsAG)^ yuiG")du, 

Js JR ^ 

We say that {3^*; e [0,r]} is an energy solution of the KPZ equation ( I2.7I I if 
there exists a constant k > such that 

E[I,AGf]<Kit-s)\\G\\l 

and 

E[{AsAG)-Al,iG)r]<^eit-s)\\G\\i 

for any < s < t < T, any e e (0, 1) and any G e 5(R). We say that a weak 
solution {yut e [0,r]} is a stationary solution if for any t e [0,r] the S'{R)- 
valued random variable yt is a white noise of variance xip)- Now we are ready 
to state the main result of this article. 

Theorem 2.4 (Equilibrium fluctuations). The sequence of proceses {{yr',t G 
[0,T]};n e N} is tight in X>([0, T], 5'(K)). Moreover, any limit point ofyi' is a 
stationary energy solution of ( 12.71 ). 

An immediate consequence of this result is the existence of weak solutions 
of the KPZ equation. Let yt be a limit point of 3^". Since the measure Vp is 
invariant under the evolution of 77", for any fixed time t e [0,T] the 5'(R)- 
valued random variable yt is a white noise of variance x{p)- The following 
corollary states some properties of the sample paths of the process yt- 

Corollary 2.5. For any limit point {yt]t £ [0,r]} of {{yr;t e [0,r]};n e N}, 
there is a finite constant c > such that the process {At;t e [0,T]} defined as 
above satisfies the moment bound 

IE[A,t(G)2] <cK-s|3/2||G||2. 

Moreover, for any 7 e (0,1/4) and any G G 5(R) the real-valued process 
{yt{G);t e [0, T]} is Holder-continuous of order 7. 

The only rigorous result about existence of solutions of ( 12.71 ) we know is the 
work of Bertini and Giacomin [7|. Let us describe this result in a precise way. 
Let C+ (K) be the set of positive, continuous functions / : M ^ R. We say that a 
process {Zt:,t e [0,T]} with trajectories in C([0, T]; C+(M)) is a mild solution of 
the stochastic heat equation 

dZt - ^^AZtdt + al3"{p)\ 4r^ZtdWt (2.10) 

if the process Zt satisfies the integral equation 

Zt = Kt*Zo- f Kt-s* ZsdW s, 
Jo 
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where Kt{x) — {2'Kip'{p)t)^^^^ exp{— x^/2(p'(p)t} is the heat kernel and * denotes 
convolution. We say that a process {ht; t e [0, T]} is a Cole-Hop f solution of ( I2.6I I 
if ht = —Lp'{p)/al5"{p) log Zt for any t e [0,T], where Zt is a mild solution of the 
stochastic equation ( I2.10I I. Defining yt = Vht in the distributional sense, we 
say that {3^t;t e [0,T]} is a Cole-Hopf solution of the KPZ equation (US if 
{ht; t € [0, T]} is a Cole-Hopf solution of ([Sjll. 

An important issue raised in [7| (see the remark after Theorem 2.1 in that 
article) is wheter a Cole-Hopf solution of ( I2.7I I actually satisfies this equation in 
any meaningful sense. Since Cole-Hopf solutions of 02. 7D arise as scaling lim- 
its of density fields in the weakly asymmetric exclusion process, Theorem 12.41 
combined with the results in El has the following consequence. 

Theorem 2.6. The Cole-Hopf solution of ( 12. 7D with initial distribution given 
by a spatial white noise of variance x{p) is an energy solution of 02.71 1. 

2.4. Current fluctuations and growing interfaces. For each x £ Z and 
any t e [0, oo), let Jp{x) be the current of particles through sites x and a; + 1 up 
to time t. That is, J"{x) counts the number of particles passing between sites 
X and X + 1. The local conservation of the number of particles plus the fact 
that the jumps are only to nearest neighbors are equivalent to the continuity 
equation 77" (x) — tiq (x) — J" (x — 1 ) — J" (x). In other words, "what we have now 
minus what we had before is equal to what came in minus what went out". The 
height function /i" = {ht{x);x e Z} is defined by 

Ki^) = \ jr(o), ^ = (2.11) 

The continuity equation shows the relation ft," (x) = h^ {x) + J" (x) . Notice that 
the process 77" corresponds to the increments of the process ft". In fact, 77" (a;) = 
h^{x — 1) — h^{x). We write this relation in the symbolic form = —Vh^. We 
can interpret ft" as the evolution of a particle system on which particles are 
created or annihilated, but on which particles never move. A jump of a particle 
from site x to site x + 1 in the exclusion process corresponds to the creation of 
a particle at site x, and a jump from a; + 1 to a; corresponds to the annihilation 
of a particle at site x. Therefore, in order to have a creation at site x, it is 
necessary to have ft"(x — 1) — ft"(a:;) = 1, ft" (a;) — ft" (a: + 1) = 0. And in order 
to have an annihilation at site x, it is necessary to have ft" (a; — 1) — ft" (a;) = 0, 
ft" (a;) — ft" (a; + 1) = 1. The creation rate, given this necessary condition, is equal 
to n'^Pnjxiht), where 7a:(ft") = C3;(— Vft"). And the annihilation rate is equal to 
n'^qnjx{h"). If a > 0, we say that ft" is agrowing interface model. It was in the 
context of growing interfaces that the KPZ equation 02.61 1 (and also 02.71 1) were 
introduced in the original work of Kardar, Parisi and Zhang [33|. 

Let us define the rescaled interface field as the process 0" = {6'"(a:);a: g M} 
by taking O'^ix/n) = n-^^'^{h^{x) - E„[ftJ'(a;)]) for a; e Z and (x) = 0J'([na;J/n) 
for X ^ n^^Z. We could also extend 0" to M by linear interpolation, but in our 
setting it will prove to be more convenient to define 0" as we do here. Notice 
that 

enx) - n-'/'{j't\Q) E„[jr(o)]) - 3^r(i(o,.]) 
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and in particular ei'{0) = n-^/^{Jl'{0) - E„[J("(0)]). Let Ho be the Heaviside 
function, that is, Ho{x) — l(o,oo)(2;)- Our first result corresponds to a limit 
theorem for the current J"(0). 

Theorem 2.7. Let yt be a limit point of y^, and denote by n! a subsequence 
such that — > yt- Then the process Jt(0) = 3^t(iJo) — 3^o(-ffo) is well defined 
and 

lim ef{0)^Jt{0) 

n' ^oo 

in the sense of convergence of finite-dimensional distributions. 

The notion of energy solutions of the equation ( I2.7D defined in Section 12.31 
can be extended in a natural way to the equation ( |2.6I >. We say that a process 
{Out & [0, T]} with trajectories in C([0, T],5'(K)) and adapted to some standard 
filtration { J"t; t e [0, T]} is a weak solution of ( I2.6I I if 

i) There exists a process {Bt] t e [0, T]} with trajectories in C([0, T],S' 
and adapted to {Fu t e [0, T]} such that for any G e 5(M), 



lim 



^{(Mi±fl^)^ _ . BAG). 



ii) For any function G e 5(M) the process 

Mt{G) = (0t,G) - (00, G) - 4^ [\es,G")ds~'^^^Bt{G) 



is a martingale of quadratic variation x{p)v'{p)t J G{x)'^dx. 

Let {6it; t e [0, T]} be a weak solution of (US. For < s < t < T, let us define 
the fields 

BsAG)=Bt{G)-Bs{G), 

«;.,(o) = /7J'"{( ^-'- + f ^'-' )' - 2^}G,.,....... 

We say that {9t; € [0,r]} is an energy solution of the KPZ equation ( I2.6I ) if 
there exists a constant k > such that 



n,G")duy <K{t-s) I G{xfdx 



E 
and 

E[{Bs,t{G) - Bl^^{G)f] < Ke{t ~ s) J G{xfdx 
for any < s < f < T, any e e (0, 1) and any G e 5(M). 

Theorem 2.8. The sequence of processes {{6*";^ G [0,T]};n e N} is tight with 
respect to the J-Skorohod topology ofV^Q, T];S'(R)). Moreover, any limit point 
{Ot;t e [0,T]} of this sequence is an energy solution of the KPZ equation i2.6i 
with initial distribution given by a two-sided Brownian motion of variance xi.p)- 
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3. Second-order Boltzmann-Gibbs principle 

Let / : f2 ^ R be a local function and let us define ip{p) — J fdvp. In ifTsll . 
the following theorem is proved. 

Proposition 3.1. Let us assume that a = 0. Then, for any function H e S{R) 
we have 

2n 



lim 



/n 



= 0. 



This result is the celebrated Boltzmann-Gibbs principle introduced by Rost. 
It roughly says that the fluctuation field associated to / is asymptotically 
equivalent to a multiple of the density fluctuation field y^. The idea is the 
following. Particles are neither created nor destroyed by the dynamics. There- 
fore, in order to equilibrate a local fluctuation of the number of particles, it 
is necessary to transport it to another region. The density of particles is the 
only locally conserved quantity of the system. Due to the ellipticity condition 
i), the process has good ergodic properties. Therefore, a fluctuation of a non 
conserved quantity will be locally equilibrated. If we look at the process in the 
right time scaling, the only fluctuation we will see will be the fluctuation of the 
density; other fluctuations being too fast to be observed in that scale. 

Notice that if = 0, Proposition 13.11 does not give a lot of information: 
it simply asserts that the fluctuation field associated to / asymptotically van- 
ishes. Let / : — > R be a local function such that ip^p) = 0. For ease of 
notation we also assume '4>{p) — 0. For G e S{M) and n e N we define the 
function V"G : Z ^ R by 



V"G = 



In other words, V"G is a discrete approximation of G'(a;/n). Let us define the 
field Ai with values in S' (R) as 

Notice that the field is analogous to the integral field appearing in Proposi- 
tion (STl] multiplied by ^/n and evaluated in G'. It turns out that the prefactor 
will make appear a non-vanishing limit for A". Another important obser- 
vation is that, unlike the case in Proposition 13.11 the field inside the integral 
does not converge to any limit as n ^ oo. In fact, its variance grows like n. 
Therefore, the convergence of Af to a well defined limit will be a purely dy- 
namical feature of the process 77". 

In order to state what we call the second-order Boltzmann-Gibbs principle, 
we need some notation. For e > we denote by eZ the set {ez; z e Z}. From now 
on, quantities like en are treated as if they were integers, meaning sometimes 
[en\ and sometimes [en] . 

Theorem 3.2 (Second-order Botzmann-Gibbs principle). For any function G e 

S(R), 



lim lim sup E„ \(aUG) - / V yi'{t,{x)f {G{x + e) - G{x))ds 



2n 

= 0. 
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This result is telling us "grosso modo" that the field is asymptotically 
equivalent to the field {y")^- Since is a distribution, the square (3^")^ 
should be defined through some type of regularization, which is exactly what 
Theorem l3.2l is saying. According to Assing |3|, the expression (3^")^ can not 
be defined as the limit in some sense of the regularizing sequence 3^"(ie(a;))^, 
at least not as a distribution. 

The proof of Theorem 13.21 is based on the multiscale analysis introduced in 
the article Il26ll . We will divide the proof into two parts. The basics elements 
used in the proof are the Kipnis-Varadhan inequality, a sharp estimate of the 
spectral gap of the generator L„ restricted to finite boxes and the equivalence 
between the grandcanonical and the canonical ensembles. In Section lSTTl we re- 
call the basic elements of the proof and in Section lS^ we explain the multiscale 
analysis. 

3.1. Elements of proof. In this section we recall Kipnis-Varadhan and spec- 
tral gap inequalities, which will allow us to estimate the variance of various 
additive functionals related to the process Af. Our aim is to establish Proposi- 
tion l3.5l and Proposition l3.71 In Section [3^ we will see that Theorem l3 .2l follows 
from Proposition 13.51 and 13.71 without further assumptions. Let / : O — > M be 
a function in L'^{vp) such that / fdvp = for any p e [0, 1]. We start recalling 
Kipnis-Varadhan inequality. 

Proposition 3.3 (DM El). For any T > 0, 



( f f{€)ds)'] <2QT\\f\Wn, 



where 



sup 

0<t<T ^ Jo 



l-i,n= sup {2{f,g)p-{g,-Lng)p] 



and (•, •)p denotes the inner product in L?[vp). 

This inequality was proved in the reversible case (a = 0) by Kipnis and 
Varadhan | 34 | and in the general case by Chang, Landim and 011a [ 14|. Kipnis 
and Varadhan also proved that this bound is sharp in the reversible case. 

This proposition is not very useful unless we have an effective way to com- 
pute the Sobolev norm Let us recall the equivalence 02. ID . This equiv- 
alence leads us to define the Dirichlet form V : L'^{vp) ^ M as 



Relation ( 12.11 1 now reads eon'^'D(f) < (/, -Lnf)p < eon'^'D(f) for any local func- 
tion / : — > M. For each x £ Z, let us define 

T^Af) = J (V.,.+i/(r?))'j/p(d77), 

sothat2?(/) ^Ex^'xl/). 

Take A C Z and / : ^ M. We say that supp(/) C A if /(r/) = /(^ whenever 
ri{x) = ^{x) for any x e A. Let us define Ta = cr{r]{x);x e A), the cr-algebra 
generated by the coordinates of i] in A. Notice that supp(/) c A if and only if / 
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is J'A-mesurable. The function Vxif) has the following convexity property: for 
AcZ and x e Z such that {x,x + 1} c A we have 

Vx{E[f\TA])<n,{f). 

Here and below, conditional expectations will always be taken with respect 
to the measure fp. The next estimate, known as the spectral gap inequality 
connects the variance of a local function / with its Dirichlet form !?(/): 

Proposition 3.4 ( ll38l [TtII ). There exists a universal constant Ao such that for 
any fc e N and any function f with supp(/) C {!,..., k}, such that J fdvp = 
for any p e [0, 1] M;e have 

//'di^p<Aofc2^ (3.1) 

x=l 

Now we explain how to use this proposition to estimate the norm 
Let / e L'^{vp) be such that / fdvp = for any p e [0, 1]. Assume that the 
support of / is contained in A = {1, . . . , fc} and let g e L'^i^p) be arbitrary. 
Define qa ^ E[g\J^A]- Then (/,g)p = {f,gA)p and 

fc-i fc-i 
{g, -Lng)p > eon^Vig) > cq ^ 1^x{g) > eo"^ ^ T^xigA}- 

x—l x—1 



Therefore, 

fe-i 

< sup {2(/, g)p - eon' ^ V,{g)}, 

S x=l 

where now the supremum is over functions g such that supp(g) C A. We 
will use the spectral gap inequality ( 13.11 ) in order to obtain a lower bound for 
Yjk ^xig)- But g does not necessarily satisfies the hypothesis of Proposition l3.4l 
Let us define g = g — E[g\ri^{0)]. Since / fdi^p — for any p e [0, 1], we have 
{fi9)p = {f-,9)p- For 1 < X < fc — 1, the transformation 77 — > j^^'^+i does not 
change the value of 77'^(0). Therefore, we also have T>x{g) ~ Now we 

can use the spectral gap inequality, since Proposition 13.41 applies for g. Using 
Proposition [3]4] we obtain that 

||/||^,„<sup{2(/,.9)p-|^(g,g)p}, 

g '^Qr: 

where the supremum is over functions g which are -measurable. This last 
supremum can be computed explicitly and it is equal to Xak'^if, /)p/eon^. Since 
the measure I'p and the Dirichlet form 'D{f) are translation invariant, the same 
estimate holds whenever the support of / is contained on an interval of size fc. 
More relevant to our purposes is that this estimate is additive in the following 
sense. Let /i, /2 : — > M be such that / fidi/p = for any p e [0, 1] and i = 1,2. 
Assume that supp(/i) C {1, . . . , fc} and supp(/2) C {fc + 1, . . . ,k + l} with fc, / e N. 



k-1 
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Then 

sup{2(/i + /2,.g)p - {g,-Lng)p} < sup{2{fi,g)p - eon'^^V^{g)} 

k+l-l 

+ sup{2(/2,g)p-eon2 ^ V^{g)}, 

9 x=k+l 

from where we conclude that 

ll/l + /2|P_i,„ < -^{k^{h,h)p + l^{f2, /2)p}. 

In other words, if the supports of /i and /2 are contained in disjoint inter- 
vals, Kipnis-Varadhan estimate is additive. Let us rewrite this observation as 
a proposition. 

Proposition 3.5. Let {Ai = {U + \, . . . + ki); i e 1} a collection of disjoint 
intervals in Z. Let fi : O — > M 6e such that supp(/i) C Ai and such that J fidvp = 
for any p e [0, 1] and any i. There exists a constant cq = co(eo) such that 

ft 



sup (/ E/»(''"H%??E^'(/-/^ 

0<t<T ^Jo'~~! ' ^ " 



iei ~ " iei 



An important final step in the proof of hydrodynamic limits is the so-called 
equivalence of ensembles. In order to state this property, we need to introduce 
some notation. For fc e N and a; e Z, we define 

k 

k I ' 

k 



In other words, ry*^ (x) is the density of particles in a box of size k at the right of 
X. Let / : 17 -> R be a local function and assume that supp(/) C {1, . . . , I}. For 
k > I and m e {0, 1, . . . , k}, we define flk = {0, 1}{^ - -M and 

fife.m = {■>] e fife; E '^(^) = "^}- 

We define the measure I'k^m as the uniform measure in Qk,m- Notice that 
Vk^m is also equal to the Bernoulli measure i^p restricted to ilk and conditioned 
to the set flk,m- We define the function 

"Hik^x) = J fdukM = E[f\ri\Q) = x] 

for any x of the form m/k, where m e {0, 1, ... , k}. The function '^{k, x) will 
play a fundamental role in the proof of Theorem 13 . 2 [ Recall also the definition 
V'(p) — J fdvp. For a local function / such that its support is not necessarily 
contained on a set of the form {1, . . . , /} there exist some positive numbers I, x 
such that supp(t^/) G {1, . . . , I}. Then we define ^'(fc, x) = E[Trcf\-q^{{)) — x]. 

Proposition 3.6 (Equivalence of ensembles). Let f : fl ^Rbea local function. 
There exists a constant Ceq = Ceq{f) such that 

l,T^/7 \ // \ 1 x(\ x) ff . .1 Ceq 
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for any a; G {0, 1/fc, . . . , 1} and any fc e N. 

This proposition is classical, and for the sake of completeness we present 
in the Appendix a simple proof for the particular case of Bernoulli product 
measures considered here. Notice that there is no mention to any particular 
density p e [0, 1] in this proposition. The following proposition explains what 
the equivalence of ensembles says for a given fixed density p e [0, 1] and for a 
function / as in Theorem l3.2l 

Proposition 3.7. Let f he a local function and let p e [0, 1] be fixed. Assume 
that ip{p) = ip'{p) — 0. Then, for any p £N there exists a constant Ceq = Ceq{f,p) 
such that 

and 

Proof The first inequality is a straight-forward consequence of Proposition [3^ 
The second one follows from Proposition 13.61 and a second-order Taylor expan- 
sion of 'ijj{x) around x — p. Let / be a local function and let p e (0, 1) be fixed. 
Then there exist bounded functions 72.* : [0,1] x [0, 1] ^ R, i = 1, 2 such that 

^{k,x) = ^(p) + -p) + ^{{x -Pf-^} 

+ n\p, x){x - pf + TZ^p, x)^^. 

□ 

3.2. The multiscale analysis. In this section we prove Theorem [3^ Roughly 
speaking, the proof has two steps. First we change the function t^/ by a func- 
tion ofT]'^{x), for a suitable k = k^. We call this step the seed. This step is 
analogous to the one-block estimate (see Section 5.4 of l.31il ). Then we define 
a sequence of scales {ki}i and we inductively change rii^^{x) by 'q^^+^{x). This 
last step is repeated until the final ki is equal to en. A careful estimation of the 
error introduced at each approximation plus a proper choice of the sequence 
{ki}i will prove the theorem. 

It will be convenient to extend a little the definition of the field AUG). Let 
us recall that we are assuming that the local function / : 17 ^ M satisfies 
'4){p) = i^'ip) = 0. For ease of notation, we assume that there is / e N such that 
supp(/) C {1, . . . , /}. Let us fix n > Z for now. Let G : Z ^ M be fixed and define 
_ff : Z — ^ R as -ffa: = n{Gx+i — Gx). We assume that 

l|G||L = -E(^-)'<+^- (3.2) 

For this function G we define A^{G) in the obvious way: 

f Y.rxf{^:)Hxds. 

Of course, our previous definition corresponds to take Hx = V"G. Let k > Ihe 
fixed. Later k will be chosen as a function of n. Our first task is to introduce 
an average of the functions Txf over a box of size k. As we mentioned before. 
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the variance of the field J2x 'Txf{ii)Hx grows Hke n. Therefore, even at this 
prehminary stage we need to use the time evolution of the system. Let us 
denote by /cZ the sublattice {kx;x e Z}. For x e fcZ, let us define 

x+k 
i—x + 1 

From Cauchy-Schwarz inequality, 

k 

{H';f<l/kY{Hx+.f 

i=l 

and we conclude that for any fc > 0, 

In what follows, we will make repeated use of this elementary inequality, with- 
out further mention to it. Let us define the rest as 

x+k 



xekZi=x+l 

= r E E nfKm-H'x)ds. 

•^0 xekZi=x+l 



(3.4) 



cSfcZ i=x+l 

Notice that for each x £ kZ the sum 

x+k 

E nf{^){H,-H^) (3.5) 

i—x + l 

has mean zero with respect to i/p for any p G [0, 1]. Since k > I , tor x,y G kZ 
such that \x — y\> k, the corresponding sums in ( 13.51 1 have supports contained 
on disjoint intervals of length at most 2k. Therefore, we can split the sum on 
the right side of ( I3.4| > into two pieces, each one of which satisfies the hypothesis 
of Proposition [5151 In this way we get the estimate 

x+k 

x£kZ •' i=x + l 
2 / J-N x+k 



14 sup (Rr'^r] < ^ E E nm{H..-H'^))\,idv) 

o<t<T n "—f J V . "-^ J 

— — x^kL i—x-^-l 

C(/)C„,fe(G)Tfc2 



where = 16co/(/, /)p is a constant that only depends on /0 and 



n 

x£kZ i=x+l 



^From now on, we do not make explicit the dependence of / of the various constants appearing. 
Therefore, constants like C(p) may also depend on /. 
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Notice that C„,fc(G) < ||G||?,„. We conclude that ^ when n — > oo as 

soon as fc^/?i — > 0. Let us define fco = 2A;, where k = k{n) is such that k^/n — > 
as n — > oo. Up to here, we have written At{G) as 

r E J2^.+JW:)H',ds (3.6) 

plus a rest that vanishes in i^(P„) as n — >^ c». Since the support of / can 
be bigger than a single point, the sums over each block of size k do not have 
disjoint support. But since k > I, if we split the sum over x e kZ into two 
alternated sums, one over x e 2fcZ and another one over x + k € 2kZ, the 
corresponding supports will be contained in disjoint intervals of size 2k. Let us 
define 

Ar = /* E E rj{v:)H',ds, 

„t x+2k 

AT = H nm)H'x+kds. 

•^'^ xG2kZi=x+k+l 

We will concentrate ourselves in Ai'"^, the computations for " being ex- 
actly the same. Notice that the support of the sum J2i '^i/('7)^^ is contained in 
the interval {x + 1, . . . ,x + 2k}. For x £ Z and k,n £N, define 

k 

^"''(^) = ^E'?"(^+*)- 

i=l 

Let us recall the definition of \E'(fc, x) and let us write 

Ar = R"r''+ f E k^{'2k-r^T^Hx))H^,ds, 



where 



xG2fc: 

x+k 



^0,n.e^ {Yl nf{lJ-)~k^{2k,^-^^'^ix))}H',. 

xe2k1 i=x+l 

The error term is defined in a similar way. Notice that k'^{2k, 7f^{x)) is 

equal to the conditional expectation of Tx+ifiv) on the corresponding box of 
size 2k. Therefore, in this integral each term of the sum over 2kZ has mean zero 
with respect to each measure Vp and we can use Proposition |3.5l to estimate the 
variance of i?" " "^. Repeating the computations done to compute the variance 
of we see that 

o<t<T frf;„' n 

A similar estimate holds for Let us define = + + 

Recall the choice ko — 2k. Putting the three estimates together and 
observing that for any x e 2kZ, + H^^^. = 2Hl^, we have proved that 

^r(G)- r V ■^{ko,rs^'^-{x))Ht'>ds + R'r{H), 



x&k 
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where 

E„|_ sup [Rf- ) \ < . (3.7) 

o<t<T n 

This decomposition is what we call the seed. What this decomposition is 
telling us, is that we can replace the weighted averages of the functions r^f 
by a function of the density of particles on a "not-too-big" block of size ka. The 
informed reader is invited to notice the parallel between this decomposition 
and the one-block estimate introduced in |[27]| . 

Now we are ready to start the multiscale argument. For e N, let us define 



Y.'^'^i'^^'is 'i^m^ds. 



We have the following estimate: 

Theorem 3.8 (Iterative bound). For any fc e N and any H bounded and of 



compact support, 

,n,2k 



sup {A7'''iG)~Al'''{G)y 

0<t<T 



^ C{p)\\G\\l„kT 



Proof. Notice that for any x e 2fcZ, 

E[^{k,r]''{x))\7j^''{x)] = E[-^{k,f]''{x + k))\r]^''{x)] = ^^{2k,r]'^'' {x)). 

This is evident from the definition of ^'(fc, ri''{0)) as the conditional expectation 
of / with respect to ry''(0). Notice as well that + H^^i, = 2H^''. Therefore, 
we can write 

2k-^{2k,Tf^{x))Hf = E[k-^{k,7^^{x))H^ 

+ km{k,r^\x + k))H^,+^\rf\x)\. 

We conclude that the functions 

F^{i-j) = k-^{k, r^^{x))H^ + k-^{k, T^^{x + k))H^^k - 2k^{2k, Tf''{x))Hf 

have mean zero with respect to each measure Vp. The support of the functions 
Fx is contained on the interval {x + 1, . . . , x + 2k}. Therefore, the functions 
satisfy the hypothesis of Proposition 13.5] By the definition of ■4}{k,ri''{x)) as a 
conditional expectation and by Proposition [321 we see that 

{Fx,Fx)p<c,,{{H^xf + {H^x+kf}. 

Since At-'^'"{G) - A'l''"{G) = JoT,xMVs)ds, using Proposition [53] we obtain 
the desired bound. □ 

For i £ N we define fc; = 2*fco- To avoid overcharged notation we write 
At'' = At'''' (G). By Theorem[3Jl we have 

G{p)\\G\\l,,TTko 



Therefore, writing 



sup [At' -At 

0<t<T 



< 



7n — 1 



i=0 
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and using Minkowski's inequality we see that 

E 
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sup (^r'™-^r'°)1 < f y]E„ 

0<t<T -I ^ 



=0 



sup (^t' -At' 

0<t<T 



1 l/2\ 2 



< 



E 

i=0 



' /c(p)||G||2„rfco2* 



< 



C{p)\\G\\l^Tko2"^ ^ C{p)\\G\\l^Tkm 



71(^2-1)2 

Notice that this last estimate only depends on km and not on fco- Fix e > 
and take m ~ log(en/fco). Recall the estimate ( I3.7D . Choosing fco — and 
putting these two estimates together, we conclude that 

E„[ sup (^J'(G)-^r''"(G))'l <C(p)||G||?„Te. (3.8) 

l-0<t<T -I 

This estimate basically finishes the proof of Theorem l3.2[ Recall the defini- 
tion ofA'^'"\Gy. 



Let us define 
We can rewrite A!t"\G) as 



en 



A'^''''{G) - 



rip) r 



•^0 KGenZ 



Notice that the constant x{p) I is not needed when iJ^^ = 0. By Proposition 
[321 we have 

Notice that 

Putting estimate ( I3.9D together with estimate ( I3.8I I and choosing = V"G, 

2n 



we conclude that 



(^r(G)-^ f {G{x + e)-G{x))y:{^,{x)fds 



< 



<G{p)\\G\\l^{Te+^}, 

which proves Theorem l3.2l Here we make two remarks about this result. First, 
we did not only obtain the convergence result stated in Theorem 13. 2[ but we 
also obtained a good control on the rate of convergence. This point will be 



22 



MILTON JARA AND PATRfCIA GONgALVES 



important in Section I4.2[ more precisely to prove that limit points of 3^" will 
be energy solutions of equation i2.7i . And second, this estimate shows that 
Theorem 13.21 holds in general for any function G : R — >^ M such that ||G||^ „ is 
finite for any n and uniformly bounded in n. 

4. Proof of the equilibrium fluctuations 

The proof of Theorem 12.41 is based on the study of some martingales asso- 
ciated to the process 3^". Let F : ^ M be a function on the domain of the 
generator L„. Dynkin's formula says that the process 

Jo 

is a martingale with respect to the natural filtration associated to the process 
ry". If F'^ also belongs to the domain of the generator L„, the quadratic varia- 
tion of m['" is given by 

{Mfn = f {L^F{ij-f ~ 2F[^-)L^F{ij:)]ds. 
Jo 

We will use this formula for F{r]f ) = 3^"(G) for G e 5(M). After some calcula- 
tions, we see that 

(4.1) 

is a martingale, where f{ri) — ^ac{r]){r]{l) — 77(0))^ and 

is a discrete approximation of G"{x/n). The quadratic variation of Ml^{G) is 
given by 

(Mr(G)) ^ f-Y. r.gniv^KKGYds, 

Jo " 

where g„(7/) = {p„r7(0)(l - + 9,i77(l)(l - r]{0))}c{rj). Looking at formula 
( 14. m . we can write ^'"(G) as the sum of four terms: the initial value yoiG), the 
martingale M^{G) and two integral terms. Starting from this decomposition, 
we will prove in Sect. 14. ll that the sequence of processes {{3^t";t > 0};n e N} 
is tight with respect to the uniform topology in 2?([0,oo),5'(R)); and then in 
Sect. l4.2l we will prove that any limit point of the sequence {{3^"; i > 0}; n e N} 
is a weak solution of the KPZ equation. 

4.1. Tightness of the density field. In this section we prove tightness of 
the sequence {{3^";t > 0};n e N}. As usual, to avoid uninteresting topology 
issues, we fix r > and we consider the processes 3^"(G) restricted to the 
interval [0, T]. We will use Mitoma's criterion ||36]| . which now we describe. Let 
X be a complete, separable metric space with metric d : X x X ^ [Q,oo) and 
let {P„;n e N} be a sequence of probability measures in T>{[Q,T],X). We say 
that {Pn\n e N} is C-tight if {P„; n e N} is tight with respect to the uniform 
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measure in V{[0, T],X). Let P be a probability measure on X>([0, T],5'(M)). For 
G e 5(R) we denote by the probability measure in ^([O, T], M) defined by 

P''iA)^P{x.{G)eA), 

where x.(G) G X'([0,T],R) is given by x.{G) = {a;t(G);t e [0,T]}. Mitoma's 
criterion says the following: 

Proposition 4.1. Let {P„;n e N} be a sequence of probability measures in 
P([0, T], 5'(M)). The sequence {P„; n G N} is C-tight if and only if{P^; n G N} is 
C-tightfor any G G S(R). 

As a consequence of this criterion, in order to prove C-tightness of the se- 
quence {{3^";< G [0,T]};n G N}, it is enough to prove C-tightness of the se- 
quence {{3^t"(G);t G [0,r]};n G N} for any function G G 5(R). It is enough to 
prove tightness for each one of the four terms appearing in H.li . First notice 
that for any G M, 

E^eMiOyriG)}] = l[E„ [exp{ — irj^{x) ~ p)G{x/n}] 

where < |G(x/n)|^. Therefore, yaiG) conveges in distribution to a normal 
random variable with mean zero and variance xip) f G{x)'^dx. We conclude 
that the sequence {yS{G); n g N} is tight (since it is convergent). The following 
proposition says that we do not need to prove C-tightness: it is enough to prove 
tightness with respect to the Ji-Skorohod topology of 2?([0, T], M). 

Proposition 4.2 (|20|). Let {P„; n G N} 6e a sequence of probability measures 
on'D{[0,T],X). assume that 

i) The sequence {P„; n G N} is tight with respect to the Ji-Skorohod topol- 
ogy of V{[0,T], X), 

ii) For any A> 0, 

lim P„( sup \x{t) - x{t-)\ >A)=0. 

n^oo o<t<T 

Then {P„; n G N} is C-tight. 

In other words, this proposition tells us that if the jumps of 3^"(G) are get- 
ting smaller and smaller with n, then Ji -tightness and C-tightness are equiva- 
lent. In our case 

sup \ynG)-yiiiG)\ < 

o<t<T 

and it is enough to prove tightness of {{3^"(G); t G [0, T]};n G N} with respect 
to the Ji-Skorohod topology in X>([0,r],M). Notice that the process Mp{G) 
also satisfies condition ii) of Proposition 14.21 The following criterion, known 
as Aldous' criterion, is very effective to prove tightness with respect to the 
Ji-Skorohod topology. 

Proposition 4.3 (Aldous' criterion [2J). Let {P„; n e N} be a sequence of prob- 
ability measures on I?([0, T],X). Let us assume that: 
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i) For any t e [0, T] and any e > there is a compact set K = K{t, e) such that 

ii) for any e > 0, 

lim limsup sup sup P„(rf(a;(T + 7, x)) > e) = 0, 

''-^O „->.oo 0<7<i5reTr 

where Tt is the set of stopping times bounded by T. Here we define x[t + 7) = 

x{t) ift + -/> T. 

Now we turn into the tightness of {M"{G)). Let us recall the definition of 
the seminorm ||G||i,„ given in Section [3^ For a function G : M — )■ M, we have 

Notice that for functions G such that G' e S{R), ||G||f „ converges to / G'{x)'^dx 
when n goes to infinity. Since the number of particles per site is at most equal 
to 1, we have the simple bound 

|(A/r(G)) - (AC(G))| < e^'\\G\\lJt-s\, (4.2) 

valid for any s, t (even random). Therefore, 

p„(|m;+^(g) - M-{G)\ >e)< 1e4(m;+^(g) - M-iG))'] 

< lE4(Af;+,(G)) - (Af;'(G))] 

This proves condition ii) of Proposition [4^3] for the martingales i\f"(G). Condi- 
tion i) follows from the fact that E„[Af"(G)^] is uniformly bounded for i e [0,r], 
71 e N. 

Tightness of the integral term 
follows from the fact that there exists a constant G = G(G, p) such that 




for any t e [0, T] and any n e N (see Il20ll or |[3T]| for more details). We are only 
left to prove tightness for 

A'iG)^ fY.Tj{r^-)V-Gds. 

It turns out that this term is the most difficult to analyze, and here we use the 
whole power of Theorem l3.2[ Notice that even to prove that A^{G) is in £^(P„) 
seems hard to prove at first glance. From estimate ( I3.8D we see that 

E„[(^r(G) < G„(p,G)ie. (4.3) 
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The constant Cn{p, G) can be chosen as equal to C{p)\\G\\\,,^. At this part of the 
argument, the only important point is that C„(/9, G) is uniformly bounded on n. 
Take = V"G and define as in ( I3.3D . Using Cauchy-Schwarz inequality 
plus the invariance of Vp under the dynamics, we obtain the estimate 

E„[y^r"(G)'] < (4.4) 
Putting these two estimates together we conclude that 

%,[Ai{GY] <G„(p,G){ie+^}. (4.5) 

Since the process 77" is stationary, the same bound holds for E„[(y^"_^^(G) — 
^"(G))^] for any s e [0,T — t]. Moreover, we can choose the constant e in a 
convenient way. However, there is a small constraint on e. At the beginning 
of the multiscale analysis in Section [3^ we took k > I. This choice, plus some 
parity considerations impose the restriction e > 21 /n. Taking e = |t — .s|^/^ we 
obtain that for any s,t e [0,T] such that |t - s| > 4P/n^ and any n e N, 

E4(^I'(G)-^^(G))'] <G(p,G)|t-s|3/2. (4_g) 

If |t — s| < AP/n'^, we can simply use Cauchy-Schwarz inequality plus the sta- 
tionarity of 77" to obtain the bound 

E„ [{AUG) - ^?(G))'] < G(p, G)\t - s\^n 

< G(p,G)2/|t-s|3/2. 

Therefore, taking a bigger constant if needed, ( 14.61 1 holds for any s, t e [0, T] 
and any n e N. Notice that there is a constant C{p) such that this estimate 
is still valid if we replace C{p,G) by G(p)||G||f „. Let us recall Kolmogorov- 
Prohorov-Centsov criterion for tightness: 

Proposition 4.4 (Kolmogorov-Prohorov-Centsov). Let {{xn{t)\t e [0,T]},n e 

N} he a sequence of continuous processes in R. Let us assume that there exist 
positive constants a, /3, K such that 

E[\Xn{t) - Xn{sT] < K\t - S\^+^ 

for any s,t e [0, T] and any n € N. Assume as well that {a;„(0); n e N} is tight. 
Then the sequence {{xn{t);t e [0, T]}, n e N} is tight. Moreover, for any 7 < a//3, 
any limit point of {{xn{t);t e [0,T]},n e N} is almost-surely Holder- continuous 
of index 7. 

Invoking this proposition, we conclude that A^{G) is tight. This finishes the 
proof of tightness of the sequence {{yt{G)-t e [0,r]},?i e N} and in conse- 
quence of {{3^t";t G [o,r]},7ie N}. 

Notice that this proposition also gives information about the limit points of 
A^{G). If {At, t e [0, T]} is such a limit point, then yl"(G) is Holder-continuous 
of index 7 for any 7 < 1/4. 

4.2. Limit points of the density field. In this section we finish the proof of 
Theorem 12.41 by showing that any limit point of 3^" is a weak solution of the 
KPZ equation. Later we also obtain some additional properties of the limiting 
points, with the aim of obtaining a uniqueness result for such solutions. 
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In Sect. I4.ll we have showed that the sequence of processes {y";n e N} is 
C-tight in ^([O, T], 5'(R)). Let {yt;t e [0, T]} be a hmit point of J^^". From now 
on and up to the end of this section we adopt the following abuse of notation: 
n will denote a subsequence for which converges to J^t. Let us recall the 
decomposition ( I4.1I I: 

yi\G) = %"(G) + I^{G) + A'liG) + Ml\G), (4.7) 

where for ease of notation we have written 

We can choose the subsequence n in such a way that the processes Af" and 
A't have limits Alt, At, well defined as processes in C([0, T], 5'(M)). By Proposi- 
tion l3.1[ the process I" converges to ^^p'ip) J^ygds. Our first task is to prove 
that Mt{G) is a continuous martingale for any G e S{R). Therefore, we need to 
prove that the martingale property is preserved by passing to the limit. This is 
not true in general under the only assumption of convergence in distribution. 
Therefore, we need an extra argument. The extra required property is uni- 
form integrability . For the reader's convenience, we recall here some elements 
of uniform integrability which will be needed. 

We say that a sequence of random variables {Xn, n e N} is uniformly inte- 
grable if 

lim sup / |X„|1(|X„| > M)dP = 0. 

The simplest criterion to verify the uniform integrability of a given sequence 
{Xn, n e N} is a moment bound: 

Proposition 4.5. If there exists p > I such that snp^E\Xn\^ < +(X), then 
{Xn] n eN} is uniformly integrable. 

And a simple criterion to show that a limit of martingales is also a martin- 
gale is the following. 

Proposition 4.6. Let {{Mt,t € [0,T]};n e N} 6e a sequence of martingales 
converging in distribution to some limit process Mt- If for any fixed time t 
the sequence {A^";n € N} is uniformly integrable, then the process Mt is a 
martingale. 

Let us go back to our problem. Since the number of particles per site is 
bounded by 1, the quadratic variation of M"(G) satisfies the deterministic 
bound (Mf(G)) < e^H\\G\\l.^. Therefore, E„[M^iG)^] = E„[(Mt"(G))] is uni- 
formly bounded, which proves that Mt (G) is a martingale. 

Thanks to Theorem l3.8[ the process At satisfies 

At{G) = \unJ^ I 3;,(z,(x))2^i^±^l— ^dxds. 

This form is slightly different to the conclusion of Theorem 13.81 where only 
a discrete summation in space is present. Considering sublattices of Z with 
different relative positions (something of the form kZ + I), we can pass from 
a discrete spatial sum to the integral above. One important consequence of 
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Theorem l3.8l is that the process yt is such that the hmit above is well defined, 
a property that does not hold for any process in I?([0, T], 5'(R)). 

Estimating the variance of (Af"(G')) using Cauchy-Schwarz inequality, we 
see that (A/"(G')) converges in L^(P„) to x{p)'p'{p)t J G'{x)'^dx. In order to 
finish the proof of Theorem 12.41 we are only left to prove that the quadratic 
variation of Mt{G) is equal to this limit. This is equivalent to prove that 
Mt(G)2 - x{p)^'{p)t\\G\\l is a martingale. We already know that Mp{G)^ - 
{M^{G)) is a martingale. Therefore, by Proposition 14.61 we only need to show 
that MJ^{Gy — (i\/"(G)) is uniformly integrable. It is immediate to show that 
the second moments of (7\/"(G)) are uniformly bounded in n, as well as the 
fourth moments of >'"(G) and I^{G). Therefore, we only need to show that 
A'f{GY is uniformly integrable. In the course of the proof of Theorem 13.81 we 
obtained a uniform bound on the norm of A" {G), which is not enough to 
prove uniform integrability. A major step of the proof was obtained using 
Kipnis-Varadhan estimate. But Kipnis-Varadhan estimate does not lead to any 
useful information about moments of order higher than 2, unless we impose re- 
strictive additional hypothesis. Therefore, we need another idea. A sequence 
bounded in is not uniformly integrable, but a sequence converging to in 

is uniformly integrable. One more time we will use the decomposition 



In ( 14. 4D we estimated the L^-norm of the last term using Cauchy-Schwarz 
inequality. Using Proposition l3.7l with p = 2, the same computations (replacing 
Cauchy-Schwarz inequality by Holder inequality) give the bound 



for any e < 1, where C{G,p) is a constant who does not depend on n, e or 
t. Therefore, A^{G) is the sum of two terms, the first one (which is .4"(G) — 
^"''"(G)) has second moment bounded by G(G, p)et and the second one (which 
is ^"'""(G)) has fourth moment bounded by G(G, p)t'^/e^ for any e < 1. From 
this observation and a good choice for e, the uniform integrability of yl"(G)^ 
follows: for t e [0, T] we have 



where G is a constant who depends only on G, p and T. Choosing e — Af^^/^ 
(actually e = for any S e (0, 1/2) would suffice) we conclude that .4"(G) is 
uniformly integrable. This proves that Mt{G) is a martingale for any G e S{R). 

Up to here, we have proved that the process {yt;t e [0,T]} is a weak sta- 
tionary solution of the KPZ equation ( 12. 7D . It only remains to prove that 
{yt-, t ^ [0, T]} is also an energy solution. This is not difficult to prove, since up- 
per moment bounds are preserved by convergence in distribution. The bound 



AUG) = --4r''"(G)} +-4r''"(G). 



E4Ar\Gf] < 



C{G,p)t^ 

£2 



¥.^[Al{Gf\(\A^{G)\ > M)] < 2E4{AnG)~Ar''{Gf] 

+ 2E„[y^^^^"(G)2l(|ytJ'(G)| > M)] 




211/2 



E[Is,t{Gf]<K{t-s)\\G\\l 
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follows from the energy estimate stated in Proposition IS.lb elow and station- 
arity. And the bound 

¥.[{As.,t{G)~AlAG)f]<C{pet\\G\\l 

follows from (|4j3]l. 

Corollarv l2.5l is follows at once from ( I4.6D and Proposition l4.4[ And Theorem 
13.21 follows by taking c = 1 and identifying the process {yt;t e [0, T]} with the 
Cole-Hopf solution of ( I2.7I I obtained by Bertini and Giacomin. 

5. Current and height fluctuations 

In this section we prove the results stated in Section 12.41 The idea is to 
exploit the formal relation n-^/^{Jl^{Q)-¥.n[Jl^{0)\) = yt'*{Ha), where the field 
y^'* is defined by 

X''*(G)=3^r(G)-:Vo"(G). 

Clearly Ho ^ and therefore this relation must be justified. We consider 
the cut-off functions G; : M — > R defined by Gi{x) — Hq{x){1 — x/1)+. This idea is 
due originally to Rost and Vares 140] and here we closely follow the exposition 
of 11291 . For each x e Z, the process J"(a;) is a compound Poisson process of 
compensator j«^+i(r/;') - jlVi,:r('?"). where 

jx,x+i{v) = n'^PnCx{v)v{x){l - v{x + 1)) 

and 

jx+i,xiv) = n'^qnCx{v)v{x + 1)(1 - 
Therefore, the process 

jn^) - f {jx.x+m:) - 3i^^M)V^ 
Jo 

is a martingale of quadratic variation 

f {rx.x+M:)+rx+iM))d'- 

Jo 

Moreover, these martingales are mutually orthogonal for different values of x. 
Using the continuity equation, we see that 

d?{0) = yl''*{Gi) + AUHo - Gi) + Ml\H^ - Gi) +Il\Ho - Gi). (5.1) 

Notice that although - Gi is not integrable, V"(i/o - G;) has bounded 
support. Therefore, the three fields A'l{Ho - Gi), Mf(iJo - G;) and ^^"(i/o - 
Gi) are well defined. A key observation here is that Gi has a small energy: 
II Gi 111,, I — l^^- According to estimates ( 14.21 1. ( I4.6I >. we see that there exists a 
constant K = K{p, T) such that 

¥.^[A^{H^-Gif]<Kl-^ 
¥.,,[M't{Ho-Gif]<Kl-\ 

We also know the dependence of K in T, but at this point the dependence 
will not be important. The following proposition, known as the energy estimate 
tells us that I"(i?o — Gi) satisfies the same estimate. 



UNIVERSALITY OF KPZ EQUATION 



29 



Proposition 5.1. There exists a constant K = K{p) such that 

M^tiGf] <Kt\\G\\ln 
for any t>0 and any G : M M such that ||G||i,„ < +00. 

This proposition follows easily from Kipnis-Varadhan inequality and corre- 
sponds to Lemma 4.2 of II38II for example. We have the following lemma. 

Lemma 5.2. There exists a constant K = K{p, T) such that 

supE„[(0r(O) < Kl-\ (5.3) 

supsupE„[(X'^*(G™)-X''*(G/))'] <Kr\ (5.4) 

m>l riGN 

The proof of ( 15.31 1 follows at once from ( 15. ID . Proposition [STT] and inequalities 
( 15.21 1 and we leave it to the reader. Estimate ( 15.41 1 follows from the identity 

y7'*{G,,,) - 3^,"'* (GO = ^r(G,„ -Gi) + Ml\G„, - Gz) +ir(G,„ -Gi). 

Proof ofTheorem \2. 71 Given Lemma |5.2[ the proof follows as in the articles 
m and mi. Let {yut e [0,T]} be a limit point of {{3^t";i e [0,T]};7i e N} 
and let P be the distribution of {yut e [0,T]}. Let i > be fixed. Recall the 
convention about the subsequence n. Approximating G; in (M) by functions 
in we see that yt{Gi) is well defined. Taking further subsequences if 

necessary, >'"(G;) converges in distribution to yt{Gi). Therefore, ( I5.4D is also 
satisfied by yt{Gi). In particular, {3^t(G;); ^ e N} is a Cauchy sequence in L^(P) 
and there exists a random variable Jt(0) such that yt{Gi) converges to Jt(0) 
when n — > 00. And from ( 15. 3D . we conclude that 0"(O) converges in distribution 
to Jt(0) as n — > 00. Convergence of finite-dimensional distributions follows in 
the same way. □ 

5.1. Convergence of height fluctuations. We start showing that 6^{x) is a 
well-defined field in 5'(M). Looking at the definition ( 12. IIP , we see that the 
field 6^ has the following representation: 

0r(x) = 0r(o)-x'(i(o,.]). 

In particular, 6'" (a;) can be written as the sum of a real-values process in- 
dependent of X and a two-sided Brownian motion. Each one of these two pro- 
cesses is well-defined in the sense of distributions. In consequence, {6^ , G) is 
well defined. 

For a function G e 5(M), let us define TqG : M ^ M as 

%G{x) = r G{y)dy. 

J —oo 

Notice that TqG e S{R) if and only if / G{x)dx — 0. Let us define as well 
A : 5(M) — > M by A(G) = / G{x)dx. The functional A is continuous with respect 
to the topology of 5(R). Define /o : R ^ M as fo{x) ^ {1 + e-^)-i. The operator 
To does not map 5(R) into itself We define T ; 5(R) ^ 5(R) as 

rG = roG- A(G)/o, 

which now maps 5(R) into itself We have introduced the operator T in order 
to relate the interface field 6*" to the density fluctuation field 3^". It will be 
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convenient to introduce the incremental interface 0"'*(.t) = 6l^{x) — Oq{x). For 
a function G e C^(R) such that / G{x)dx — 0, it is easy to see that 

{e^,G)^ynTG). 

In terms of the field 6*" *, this relation reads 

(0,"'*,G) =X"'*(rG). (5.5) 
For arbitrary G e 5(R), the following lemma shows howto compute (0" *, G). 
Lemma 5.3. There exists a process {X '*(/o); t ^ [0, T]} such that 

{er,G) = y^iTG) + A{G)y:-*{fo) 

for any G G S{R) and any t e [0, T]. 

Proof. Our first task is to prove the existence of the process *(/o)- The idea 
is to introduce a cut-off like the one introduced by Rost and Vares. But this 
time we need to do it in a smooth way. Let us consider the bump function 

C : R -> [0,oo) defined by C(a;) = cexp{-l/x(l - x)} for x e [0, 1] and C(a;) = 
otherwise. The constant c is chosen in such a way that / g{x)dx = 1. Starting 
from this bump function we define the interpolating function g : K ^ [0, 1] as 

{1, ifa;<0 
l-S;av)dv, ifO<.T<l 
0, ifl<a;. 

Notice that g'{x) — ~C{x). We define then the cut-off functions /q : M — > R as 

/'(a;) = fa{x)g{x/l). We have that e 5(R) and 

^(X) - ax) = f',{x){9{x/l) - 1) - y/o(x)C(x/0. 

Both terms in this last sum converge to in L^{Z). Therefore, ||/q - /o||i,ri 
as ^ — > oo for any n. A more careful computation reveals that the convergence 
is uniform in n. It can also be checked that A"/q converges to A"/o in L^(Z) 
and the convergence is uniform in n. From 04. ID we know that 

In SectioniUwe obtained the estimate E„[Aff(G)2] < G(p)i||G||? „, vahd for 
any function G e S{M). In particular, the sequence {M"(/o); Z G N} is a Cauchy 
sequence in i^(P„). Therefore, for any t e [0,r] there exists a limit process 
Mf(/o) such that 

lim E„[(Afr(/^)-Mr(/o)f]=0. 

By Doob's inequality, the limit is a martingale and we also have convergence 
at the process level. The same is also true for the processes I" (/q) and A^{fl)■. 

E„[(ir(/^)-xr(/n)'] < ^l^^(A^/^-A^/™)^ 
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Therefore, the hmiting processes Itifo) and .4"(/o) are also well defined. 
The convergence at the process level follows from tightness arguments, like in 
Section mi 

At this point we just define 3^"'*(/o) using ( I4.1I >: 

yrifo) = Mr(/o) +xr(/o) + AUfo). 

For G e 5(R) such that / G{x)dx = 0, the lemma follows from ( I5.5I I and 
an approximation procedure as above. By linearity, we only need to prove the 
lemma for a single function G such that / G{x)dx ^ 0. For arbitrary func- 
tions, the lemma follows from linearity. But we have already done it in the 
previous section. In fact, at a formal level, what we proved was the iden- 
tity (0"'*,(5o) = yt'*{H(i), where 5q is the 5 of Dirac at a; = (notice that 
6'"'*(0) = 0"(O)). A similar argument proves the same identity for the bump 
function C,. Notice that we are not able to show convergence at the level of 
processes for the approximations when / G{x)dx ^ 0, but only in the sense 
of finite-dimensional distributions. But this is enough to our purposes, since 
we have already proved that >'"'*(/o) is well defined as a process, and finite- 
dimensional distributions characterize processes in 2?(5'(M), [0, T]). □ 

Now we have the elements to prove Theorem 12.81 From 04. ID and Lemma 
15. 3[ we have 

(0,"'*, G) = irCToG) + ^rCToG) + Af,"(roG). (5.6) 

By Proposition I4.2[ in order to prove tightness of e [0,r]};7i e N} 

it is enough to prove tightness of the projections {(0"'*,G);i S [0,r]};7i e N}. 
But this follows as in Section l4rTl since the estimates for the processes at the 
right-hand side of l5.6l onlv depend on the regularity of -^%G — G. We have 
already noticed that {0^{x)] a; e M} converges to a two-sided Brownian motion. 
Therefore, tightness of {{6'";< e [0,r]};n e N} follows. We are only left with 
the problem of identifying the limit points of 0". As in Section |4l2l we denote 
by n a subsequence such that all the processes above converge in distribution 
to the corresponding limit. Notice that 1\'^%G = %KG for any G e 5(M). 
Therefore, 

' y^{^^,%G)ds= f y^{ToA^G)ds, 
Jo 

and J"(7oG) converges to ^(p'{p) Jq{&1\ G")ds in distribution as n ^ oo. Since 
V^ToG = ToV^'G, the process MpiToG) converges to Mt{G), where Mt is a 
space-time white noise of variance x{p)^'{p)- It remains to identify the limit of 
in terms of 9t. It is enough to write ^" '^"(G) in terms of 6'" . We have 







9^{x)-e^{x + e) 



e 

and therefore the limit process At satisfies 



A(roG) = lini ^-^{ r^^^^^^-'^^^^ Y^^}Giy)dydxds. 
<^^oJo JrJx 2e l\ e J e ) 

Notice the Wick renormalization factor x{p)/^ in this formula. This factor is 
now needed because / G{x)dx ^ 0. 
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Passing to the hmit in ( I5.6I I, we see that any hmit point {0t;t e [0, T]} of the 
sequence {6'"; i e [0, T]}; n e N} satisfies 

{Ot, G) = {00, G) + ^ / {es,G")ds + + Xt(G), 



where 64(0) = At{%G) and A^t(G) is a martingale of quadratic variation 
x{p)^'{p) J G{xYdx, which proves that {0f ; t e [0, T]} is a weak solution of ([2^. 
In a similar way, we can prove that {9t; t e [0, T]} is an energy solution of ( I2.6D . 

Appendix A. Equivalence of ensembles 

In this appendix we prove Proposition [3]6] Let / : O — > K be a local function 
and assume without loss of generality that there exists I e N with supp(/) c 
{1, . . . , /}. Our first observation is that / is a linear, finite combination of func- 
tions of the form n^eA vi^)^ where A c {1, . . . ,1}. Since the thesis of Propo- 
sition 13.61 and 13.71 are preserved under linear transformations, it is enough to 
prove the propositions for functions of the form Hi eA vi^)- We call these func- 
tions monomials. The random variable ?7'^(0) and also the measures Vp, Vk,m 
are exchangeable, in the sense that they remain unchanged under a permuta- 
tion of the random variables {??(!), • ■ ■ , ?7(fc)}. Therefore, it is enough to prove 
the propositions for functions of the form f{ri) = 77(1) • • • ri{l) with I e N. Fix 
I e N and take k > I. Let us recall the definition ^'(fc, x) ~ E\j\rf(^) — x] for x 
of the form m/k, with m = 0, 1, . . . , fc. The conditional expectation E[f\ri'''{0)] is 
easy to compute: 

Let us call ^ the first product in the last display. Notice that Uk^i is uniformly 
bounded in k, and it converges to 1 as fc — > 00. Developing the second product 
we have the expansion 

/ 

^{k,x)^auJ2^^'~'- (A.2) 

i=0 

The coefficients Pi do not depend on m or fc. Therefore, all the powers of order 
smaller than n — 1 in iA.2i are at most of order 1 /fc^, uniformly in x (recall that 
< a; < 1). Therefore, there exists a constant Ki such that 



sup |*(fc,x) - afc,z(poa;' + — M|< 



k ^1 - fc2 ■ 



Now we just need to compute pa and pi. The constant pa is equal to 1, since 
in the last product in dA.ll l each factor is monic. By the same reason, the 
coefficient pi/k is equal to minus the sum of the roots of the aforementioned 
polynomial. Therefore, pi — - l)/2. Up to here we have proved that 

sup |*(fc, x) - ak,i (x^ ~ ^^^2;"-!) I < (A.3) 

Now we turn into ak,i. It is easier to expand l/ak.i- In fact, we have that 
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wherer(fc) is bounded in fc. Using the expansion (1 — e + 0(e^)) ^ — l + e + 0(e^), 
we see that 

for another function f(fc) bounded in k. Putting this asymptotic expansion for 
flfc ; back into dA.Bl l, we conclude that there is a constant Ceq which only depends 
on I such that 



sup 



a;) 



1 



2fc 



2fc 



< 



fc2 



Remember that for this particular choice of / we have ip{x) = and (p"{x) — 
l{l—l)x''~'^. Replacing above by i^(a;) and /(Z— by xip"{x), Proposition l3.6l 
is proved. 
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